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ABSTRACT 

Many exoplanets have been observed close to their parent stars with orbital 
periods of a few days. As for the major satellites of the Jovian planets, the 
figure of these planets is expected to be strongly shaped by tidal forces. How- 
ever, contrarily to Solar System satellites, exoplanets may present high values for 
the obliquity and eccentricity due to planetary perturbations, and may also be 
captured in spin-orbit resonances different from the synchronous one. Here we 
give a general formulation of the equilibrium figure of those bodies, that makes 
no particular assumption on the spin and/or orbital configurations. The gravity 
field coefficients computed here are well suited for describing the figure evolution 
of a body whose spin and orbit undergo substantial variations in time. 

Subject headings: celestial mechanics — planetary systems — planet-star interactions 
— planets and satellites: general — planets and satellites: physical evolution 



-3- 

1. Introduction 

All the main satellites of the giant planets in the Solar System are in synchronous 
rotation, and locked in a nearly zero-obliquity Cassini state; their orbits are nearly circular, 
with orbital periods less than 16 days, and lie in the equatorial plane of the central planet^. 
These configurations are the result of the long-term evolution of the spin and orbits due 
to tidal forces raised by the central planet (e.g. Correia 2009). The influence of the tidal 
deformation on the shape of the satellites is also appreciable, because the satellites' figures 
are supposed to be in hydrostatic equilibrium with the tidal and centrifugal potentials (e.g. 
Schubert et al. 2004). This assumption results from the fact that their interiors are either 
hot (lo and Europa) or plastic (icy satellites, e.g., Ganymede, Callisto, and Titan). 

About 40% of the presently known exoplanets are found in orbits with periods below 
16 days, and at least more than 25% of them are not alone in their systems^. As for the 
satellites above, we can assume that these planets acquire an equilibrium shape dictated by 
its internal gravity and a perturbing potential. However, due to planetary perturbations 
from additional companions, close-in exoplanets may present non-zero obliquity (e.g. 
Laskar & Robutel 1993; Levrard et al. 2007) or eccentric orbits (e.g. Correia & Laskar 2004; 
Mardling 2007). In addition, they may also be captured in spin-orbit resonances different 
from the synchronous one (e.g. Goldreich & Peale 1966; Rodriguez et al. 2012). 

Previous works have shown that non- synchronous satellites present a different figure 
(Zharkov & Leont'ev 1989; Giampieri 2004). In addition, if the orbit has some eccentricity, 
the low order gravity coefficients can be separated into static and periodic components 
(involving different time-scales), and only the static part contributes to the permanent 



^http://ssd.jpl.nasa.gov/ 
^http : / /exoplanet . eu/ 
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deforination (Rappaport et al. 1997; Giampieri 2004). Therefore, since exoplanets may 
evolve in very eccentric orbits, their figures are most hkely different from what we observe 
in Solar System satellites. 

In this Letter we generalize previous studies to any eccentricity value. We also include 
the effect of an arbitrary obliquity, that may change completely the gravity field coefficients. 
We adopt a vectorial formalism, which is different from the traditional expansion of the 
gravitational potential in spherical harmonics. 



2. The gravitational potential 

The gravitational potential of a body of mass ?7i at a generic point P is given by: 



r — r 



(1) 



where G is the gravitational constant, r is the position of P, and r' the position of a mass 
element dm with respect to center of mass of the body. Assuming r' <^ r, we can develop 
the potential limited to the second order, which gives (e.g. Tisserand 1891; Smart 1953): 



V{r) 



Gm 3G 



2^3 



f ■ X ■ r 



-tr(X) 



(2) 



where r = r/r. 



III 1 12 Iiz 

Ii2 I22 -^23 l"^/ 

-^13 -^23 -^33 

is the itertia tensor, and tr(X) = In + I22 + -^33 its trace. We chose as reference the 
non-inertial frame (i,j,k) fixed to the planet, where k is the direction of the spin axis 
(Fig.l). 
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Fig. 1. — Reference frames for the definition of the reference angles, n is the normal vector 
to the orbital plane of the perturbing body at r^. The reference frame (i, j, k) is fixed with 
respect to the planet's figure, k being the rotation axis, that is tilted by an angle e (the 
obliquity) with respect to n. The rotation angle, 9, and the position of the perturber, ip, 
are measured from the line of nodes between the equatorial and orbital planes, along their 
respective planes. Vt is the rotation rate. 
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The above potential is also often expressed alternatively using the gravity field 
coefficients J2, C2n, S2n, as 



V{r) 



Gm GmR 



2 






C2oA(k-r) 
+ C2i3(i-f)(k-f) + ^2i3(j-f)(k-r) 
+ C22 3((i-fr-(j-f)^) + 522 6(i-r)(j-f) , (4) 



where ^2(2;) = (3x^ — l)/2, and 



-^33 ~ 2 (-^11 "'" -^22 J 



C21 = 5^ , 5*21 = 5^ , (6) 

_ I22- hi r. _ -^12 /«N 



3. The perturbing potential 

The mass distribution in the body, characterized by the inertia tensor (Eq. 3), is a 
result of the self gravity, but also of the body's response to any perturbing potential. 

For a planet with rotation rate f2 = fi k, a mass element dm is also subject to the 
centrifugal potential (e.g. Goldstein 1950) 

K(r') = -^in X r')^ = -In^r'^ (l - (k ■ r')^) . (8) 

In addition, since the planet moves around a central star of mass m^, we also need to 
consider the tidal potential (e.g. Lambeck 1980) 

^*(r',rJ = -^^f-yP2(r'-r.), (9) 
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where r^, is the position of the star with respect to the planet center of mass. In the frame 
of the planet, we can express r^, as (e.g. Correia 2006): 

cos e sin ip sin 9 + cos ip cos 9 

cos e sin t/' cos 6' — cos '?/' sin ^ , (10) 

— sin e sin ifj 

where e is the obliquity of the planet (the angle between the spin axis and the normal to 
the orbit), 9 is the rotation angle, ip = u + v i^ the angle between the line of nodes and the 
direction of the star, u is the argument of the periapse, and v is the true anomaly (Fig. 1). 

The resulting perturbing potential is then given by 



Vp(r\v,) = V,{v') + Vt{r\v,) 



;ii) 



that can also be rearranged as 



V,{v\v,) 



where 



with 



and 
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(12) 



(13) 



(14) 



(15) 



The last term in the above equation is directly obtained from expression (10), where ^ 
denotes the transpose. 



4. Equilibrium figure 

The perturbing potential (Eq. 11) acting on the planet, deforms it, and modifies the 
external gravitational potential (Eq. 1). A most convenient way of defining this deformation 
is through the Love number approach (e.g. Love 1927; Peltier 1974), in which the body's 
response is evaluated in the frequency domain. As long as the distortions are small, we 
can simplify the problem by ignoring the small interaction terms between the centrifugal 
and tidal potentials (Zharkov & Trubitsyn 1978). The additional gravitational potential, 
A\^(r), due to the deformation of the planet in response to the perturbing potential, is then 
given at the planet's surface (r = R) by (e.g. Lambeck 1980): 

AViR) = k2VpiR,r,) , (16) 

where k2 is the second Love number for potential. In general, k2 is a complex number, 
which depends on the frequency a of the perturbation. \k2\ gives the amplitude of the tide, 
while the imaginary part gives the phase lag between the perturbation and the maximal 
deformation. 

If we assume that the planet behaves like a Maxwell body^ with homogeneous density 
p we have (Henning et al. 2009): 

k2{a) = kf\±^ , (17) 

1 + 1 CTTb 

where kj is the fiuid second Love number, r^ = v/fi is the relaxation time, and 

U f ^ ^ 19/i-R ^ kf 



1 -4- 1 = ns 

Ta V 2Gmp) k2{oo) 

kf is constant for a given body and corresponds to its maximal deformation (for a static 

perturbation a = 0, and hence k2 = kf). For an homogeneous sphere kf = 3/2, but more 



^A Maxwell body behaves like an elastic body over short time scales, but fiows like a fiuid 
over long periods of time. It is characterized by a homogenous rigidity fi and viscosity v. 
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generally kf can can be obtained from the Darwin- Radau equation (e.g. Jeffreys 1976): 



733 2 2 i-k, \ 

^ = 3 (' - 5 VTTfc;j ■ ''«> 

In Figure 2 we plot the dependence of \k2\ with the tidal frequency for the Earth. Since 
static and elastic stresses involve very different time scales, there is no conflict between the 
two types of responses (Fig. 2). We can thus assume a static response for the centrifuge 
potential (Eq. 8) and an elastic one for the tidal potential (Eq. 9), which depends on the 
varying position of the star with respect to a point at the planet's surface. 

Combining expressions (2, 12, and 16) we thus have 

I = lo + k2lp{R) , (20) 

where Xq oc I corresponds to the moment of inertia with spherical symmetry, I being 
the identity matrix. It becomes then straightforward to compute the gravity coefficients 
(Eqs. 5—7), for which Xq is irrelevant. 

For the centrifuge contribution, all terms in Xc are zero except 133 (Eq. 14), thus the 
only non-zero gravity coefficient is J2: 

.j.,,(o,m^.,|^. pi) 

For the tidal contribution, the inertia matrix is permanently modified (Eq. 15), meaning 
that all gravity coefficients are non-zero. The global contribution to the gravity coefficients 
is given by the sum of the centrifuge and tidal contribution (Eq. 13), i.e., 

J2 = 4 + 4 , (22) 

with 



, _ A^m^ f R^ ^ 
2 m \ r^ 



3 

1 — - sin^ £^ (1 — cos 2^) 



(23) 
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Fig. 2. — Dependence of the amplitude \k2\ with the tidal frequency a for the Earth. We 
use a visco-elastic model (Eq. 17) with v ~ 10^^kgm"^s~^ (Cizkova et al. 2012), k2 = 0.299, 
and kf = 0.933 (Yoder 1995). 
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while for the remaining coefficients, the tidal contribution is the total one: 

3 



c. 



21 



2 m 

2^ ■ / 
— cos - sm 

2 ^ 



sine 



cos e sin 6 

e 



2ip) + sin^ - sin(^ + 2^;) 



(24) 



'21 



2 m 



sine 



cos e cos 9 



cos" - cos(6' — 2ip) 



sm - cos 



{d + 2^) 



(25) 



a 



22 



8 m \r^ 



sin e cos 26* 



+ 2 cos^ - cos 2{e -^) + 2 sin^ - cos 2(^ + ?^) 



(26) 



S 



22 



8 m \r^ 



sin e sin 2^ 



+ 2 cos^ - sin 2(^ -^) + 2 sin^ - sin 2(^ + ?/;) 



(27) 



Note that, unless the planet evolves in a circular orbit, we cannot replace k2 by kf in 
expression (23) for the terms independent of ip, because the distance to the star also varies 
with ip as (e.g. Murray & Dermott 1999): r^ = a(l — e^)/(l + ecosi/), where a and e are 
the semi-major axis and the eccentricity of the orbit, respectively. 

5. Permanent deformation 



In the previous section we assumed for the tidal perturbation an elastic response, but 
this is not completely true for all harmonics. Indeed, although the position of the star with 
respect to the planet surface may not be constant, when we average its motion over one 
orbital period, some perturbations do not average to zero, and the planet can assume a 
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different permanent figure. In order to identify tlie static liarmonics, we need to perform an 
expansion of the true anomaly v in series of the eccentricity e and mean anomaly M: 



^ = :^i:G(.,^)e*''-*', (28) 



(J= — oo 

and 

e 



i(ke-2i/) -|^ +°° 



rt a^ 



Y, i^(g,e)e'('^^-2^^) , (29) 



where the functions G{q,e) and H{q,e) are power series in e (Tab. 1), and k and 2q are 
integers^. Static perturbations thus correspond to frequencies for which a = k6 — 2qM = 0. 
They can occur whenever /c = g = 0, or if there is a commensurabihty between the rotation 
rate and the orbital motion (spin-orbit resonances). 



5.1. Fast rotating planets 

For an arbitrary rotation rate (such as the rotation of the Earth), only terms with 
k = q = will contribute to a permanent deformation of the planet. By replacing equations 
(28) and (29) in the expressions of the gravity coefficients (Eqs. 23—27), and averaging over 
the rotation angle 6 and the mean anomaly M, we get that all gravity coefficients become 
zero, except J2: 

(30) 



1 3.2 

1 sm e 

2 



< J2>= if — - G{0,e 
2 m \a J 

where G(0,e) = (1 — e^)~^^^. Although for eccentric orbits the distance to the star 
constantly varies with ip, there is a permanent perturbation along the direction of the two 
bodies whose average is not zero. 



^We have retained the use of semi-integers for q for a better comparison with previous 
works. 
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Table 1: Hansen coefficients G{q,e) and H{q,e) to e^. 



Q 


G(g,e) 


Hiq,e) 


-1/1 


fe^ + |e^ 


24^ 


-1/2 




48 '^ 


0/1 


1 + fe^ + f e^ 





1/2 


2e -r ^qC 


2^ -r igc 


1/1 


9 2 , 7 4 
4^ -r 4e 


1 _ 5 2 , 13e4 


3/2 


53 „3 
16 '^ 


7„ 123 3 
2*^ 16 ^ 


2/1 


77 4 
16 '^ 


17^2 115„4 
2 '^ 6 '^ 


5/2 




845^3 
48 ^ 


3/1 




533 „4 
16 ^ 



The exact expression of these coefficients is given by (^(g, e) = ^ /J^ (^) exp(i2qM) dM and 
H{q,e) = 1/; (f)'exp(i2i.)exp(i2qM)dM. 
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5.2. Spin-orbit resonances 

When the spin is tidally evolved, it may be captured in a spin-orbit resonance, for 
which Q = pn (e.g. Goldreich & Peale 1966; Correia & Laskar 2009). The most common 
outcome is the synchronous resonance {p = 1), observed in the Solar System for the 
majority of the main satelhtes. However, non- synchronous configurations are also possible, 
as it is the case of Mercury (p = 3/2) (Colombo 1965). 

For synchronous rotation, the planet acquires a permanent deformation along 
the direction that always point to the perturber (e.g. Ferraz-Mello et al. 2008). For 
non- synchronous resonances, in average, one direction always points to the perturber at 
the periapse (Goldreich & Peale 1966). Thus, bodies with some rigidity can also acquire a 
permanent deformation along this direction. 

In order to obtain the contribution to the gravity coefficients (Eqs. 23—27), we average 
again over 6 and M. However, since the rotation rate is now in resonance, we have that 
(f) = 9 — pM is constant, and therefore we must retain the terms with argument (g = pk/2) 
in the expansions (28) and (29). The contribution to Jg is still given by expression (30), 
but the remaining gravity coefficients become: 

< C21 > 



< ^21 > 



kfm^ fRV 
— - — — sine 
2 m \a J 


"g(p/2. 


H{p/2, e) cos^ ^ sin(0 - 2u) 


H{-p/2) sin^ - sin(0 + 2uj) 


kfm^ /rY 
' ' sine 
2 m \a / 


'g(p/2. 


H{p/2, e) cos^ - cos((; 


h-2uj) 



(31) 



2^ 



H{-p/2, e) sin^ - cos(0 + 2uj) , (32) 
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<C,,.> = -^— - 



kf m^ f R 
m \ a 



G{p, e) sin^5 cos 20 



2H{p, e) cos^ - cos 2(0 — u) 
2H{—p, e) sin^ - cos 2(0 + cu) 



(33) 



< Soo > 



kf m^, / R 
8 m \ a 



G{p, e) sin^e sin 20 



+ 2H(p, e) cos^ - sin 2(0 — cu) 
+ 2H{-p,e)sin^^sm2{(f) + uj) 



(34) 



Notice that the coefficients < C21 > and < S21 > can only be different from zero for 
"integer" spin-orbit resonances (1/1, 2/1, 3/1, etc.), since the Hansen functions G{p/2,e) 
and H{p/2, e) are not defined when p is an half-integer (Tab. 1). 

When the argument of the periapse, u, is also a fast varying angle {io ^ Tq7^), which 
is often the case for solid close-in planets and satellites, the resonant angle becomes 
7fc = — ku, with k = 0, ±1, or ±2 (Correia & Laskar 2010). We can therefore also average 
the gravity coefficients over u, only retaining terms in 7^ for a given k value. For moderate 
obliquity, the dominating term is 71 = — w. Thus, we get 



< C21 >=< 521 >= 



< Coo > 



kfTTl^ /RV rjr . i£ ^ 

— H{p,e) COS -cos 271 , 

4 m \a J 2 



(35) 



(36) 



< ^22 > 



— J- — - f — ) if (p, e) cos^ - sin27i 



4 m 



e 
2 



(37) 



For damped librations 71 is constant. Addopting 71 = 0, i.e., the projection of the i-axis in 
the orbital plane points to the star at periapse, we further get cos 271 = 1 and < 5*22 >= 0. 
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For a planet with the spin-axis normal to the orbit [e = 0), and truncating the series 
H{p,e) to e^, we retrieve the same results as in Giampieri (2004). In Table 2 we compute 
the gravity coefficients for the terrestrial planets and main satellites in the Solar System 
and compare with the observed values. We obtain a good agreement in all situations except 
in the cases of Mercury, Venus and the Moon. The observed values in these three situations 
correspond to fossilized values acquired when their rotations were much faster than today 
(Touma & Wisdom 1994; Correia & Laskar 2001, 2012). 

6. Conclusion 

Using a vectorial formalism, we derived the gravity field coefficients of a planet 
in hydrostatic equilibrium with the tidal and centrifugal potentials. We have made no 
particular assumption on the inertia tensor, so our results are valid for any rotation rate, 
obliquity and orbital eccentricity. In particular, they allow us to compute the shape of the 
planet for an arbitrary spin-orbit resonance. 

Combining expressions (21), (30) and (36) we get 

C22 3H{p,e) cos\e/2) 



{31 



J2 (2p)2+(6-9sin2e)(l-e2)-3/2 ' 

which may considerably differ from the ratio 3/10 observed for the satellites in the Solar 
System (p = 1, e = 0, e = 0) (Table 2). Therefore, if we are able to measure the ratio C22/ J2 
for exoplanets in eccentric orbits, for instance, by detecting differences in the light curve 
at each transit (e.g. Barnes & Fortney 2003; Ragozzine & Wolf 2009), we can determine 
in which spin-orbit the planet is locked, assuming zero obliquity. Inversely, if we assume 
synchronous rotation, we can infer the obliquity of the planet. 

Our work is also important for future studies on the long-term evolution of planets 
and satellites. Indeed, as the shape of the planet changes from one spin-orbit resonance to 
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Table 2: Low order gravity field coefficients (xlO 



Body 


kf 


< J2 > 


J2 


< C22 > 


C22 


Mercury " 


0.928 


0.534 


50.3±0.2 


0.068 


8.1±0.1 


Venus ^ 


0.928 


0.0516 


4.46±0.026 


- 


0.539±0.008 


Earth ^ 


0.933 


1063.5 


1082.6 


- 


1.57 


Mars^ 


1.20 


1830.3 


1960.5±0.2 


- 


-54.73±0.02 


Moon^ 


1.44 


8.968 


203.8 ±0.6 


2.675 


22.4±0.1 


10^= 


1.29 


1834.5 


1859.5±2.7 


546.8 


558.8±0.8 


Europa '^ 


1.04 


431.0 


435.5±8.2 


129.5 


131.5±2.5 


Ganymede '^ 


0.80 


125.97 


127.5±2.9 


37.92 


38.26±0.87 


Callisto^''^ 


1.11 


34.3 


32.7±0.8 


10.3 


10.2±0.3 


Rhea'^ 


1.25 


780.8 


794.7±89.2 


234.9 


235.3±4.76 


Titan ^ 


1.00 


32.9 


33.5±0.6 


9.86 


10.02±0.07 



The comparison is done between the values obtained through the averaged equations 

(Eqs. 21, 30, 36) and the values given by the observations. < C22 > is computed only when 

the rotation is trapped in a spin-orbit resonance. References: " Smith et al. (2012); ^ Yoder 

(1995); " Schubert et al. (2004); '^ Anderson et al. (2001); " less et al. (2007); f less et al. 

(2010). 
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another, the capture probabilities are considerably modified. As an example, if one suppose 
that the Moon acquired its present figure in the past when it was closer to the Earth, we 
conclude that the Moon was not synchronous at the time, since C22/ J2 = 0.11 (Tab. 2). If 
the rotation was in a 3/2 spin-orbit resonance (Garrick-Bethell et al. 2006), using equation 
(38) we can determine that e ~ 0.17, and from expression (36) that a ~ 25 i?©. 

We acknowledge support from FCT-Portugal (PTDC/CTE-AST/098528/2008, 
SFRH/BSAB/1148/2011, PEst-C/CTM/LA0025/2011), and FAPESP (2009/16900-5, 
2012/13731-0). 
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